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Abstract 

We solve the Faddeev equations for 11 Li (n+n+ 9 Li) using hyperspher- 
ical coordinates and analytical expressions for distances much larger than 
the effective ranges of the interactions. The lowest resonances are found 
at 0.65 MeV (| + ,| + ,| + ) and °- 89 MeV (| ± ) with widths of about 
0.35 MeV. A number of higher-lying broader resonances are also obtained 
and related to the Efimov effect. The dipole strength function and the 
Coulomb dissociation cross section are also calculated. 
PACS numbers: 21.45. +v, 11.80. Jy, 21.60.Gx 

Introduction. The new class of nuclear states, called halos, are presently 
intensively investigated pi. Their structure and behavior in reactions differ 
substantially from those of ordinary nuclei. They are unusually large and loosely 
bound. Most of them are presumably not identified yet. Large electromagnetic 
dissociation cross section, or a concentration of 1 — strength function at low 
energies exemplify the unusual character of these nuclei. Descriptions as two- 
and three-body systems have been successful ^ |3) . The continuum structure 
of halo nuclei is still almost unknown even for the most studied halo nucleus 
n Li. However almost inevitably low-lying resonances must be present, since the 
Borromean nature implies that the two-body subsystems are close to threshold 
and perhaps nearly fulfill the Efimov condition j|, |j| . 

A recent experiment, where protons are scattered on 11 Li, was interpreted 
as evidence for a l - -resonance in 11 Li at 1.0 ± 0.1 MeV above the two-neutron 
threshold with a width of about 0.75 ±0.6 MeV M. An other interpretation 
of the same experiment is given in terms of the shake-off process M . Reaction 
experiments also indicate an excited 11 Li-state at about the same energy || ||. 
The measured §, |, [n| and computed O, |lj, [l3| 1 -strength functions are 
available in the literature, where the absence of calculated 1~ resonances is 
quoted without details. 

Clearly reliable three-body continuum calculations are desirable for n Li, but 
this is not easy due to both technical difficulties and lack of information of 10 Li. 
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However, recently a method dealing with the necessary large distances was for- 
mulated for s- waves || Q and developed and applied for three a-particles and 
two neutrons surrounding a spin zero core |l5| . The purpose of the present let- 
ter is to extend the method to finite core-spin and provide continuum structure 
calculations for 11 Li. 

Method. The method is described for bound states in ^, [|, [l4) and for con- 
tinuum states in [jl5). We shall here briefly sketch the necessary generalizations. 
The k'th particle has mass rrik, charge eZk, spin Sfc and coordinate r^. The 
two-body potentials are Vij . The spins Sj and Sk are coupled to s which in turn 
coupled to Si gives the total spin. The corresponding wave function is Xs ■ The 
total wavefunction is expanded in a complete set of hyperangular functions 

' n— 1 v ' 

where each of the three components 0„ is expressed in the corresponding set of 
hyperspherical coordinates (p, fli) — (p, a>i, fl X i, £l y i). They satisfy the angular 
part of the Faddeev equations with the eigenvalues X n (p). 

For Borromean systems we select those solutions \fv to eq.(|l|) where the 
large-distance (p — » oo) boundary conditions for /„ ' are given by 

fi n 'Hp) "> ^"'M - 5«,n'4 +) M - ( 2 ) 

where k 2 = 2mE/h 2 and Fn are related to the Hankel functions of integer 
order 

fWM,^<» 2 M. ( 3 ) 

where K n {K n + 4) = X„(p — oo). The continuum wave functions are orthogonal 
and normalized to delta functions in energy. 

By diagonalization of the S'-matrix we obtain eigenfunctions and eigen- 
phases. Resonance energies and widths can be obtained by use of the complex 
energy method, where the Faddeev equations are solved for E = E r — iT/2 with 

the boundary condition /4" ' oc 8 n . n i ^J~^^ H^ +2 (np) These solutions corre- 
spond to poles of the S'-matrix 

Large-distance behavior. We shall now specialize to the system of interest, 
i.e. two neutrons [k = 2,3) and a core (k = 3) with spin s c . The three 

(i) 

independent orbital s-wave components, s ,,i = 1,2,3, are characterized by 
Ixi = 0, l yi = L and the spin wave functions x£Lo> X^L 8c _i/2! ^f 3 L Sc +i/2- 
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All other partial waves decouple for large p over a distance scale defined by 
the range of the interactions. The s-waves, however, are coupled and feel the 
interactions over a distance defined by the scattering lengths. The remaining 

(2) (3) 

two s-wave components are determined by antisymmetry, i.e. (j>x, s = ~ < Pl S3 i 
f>L 3 = 4>l\ 2 - The l as ^ combination, 4>^\=ii is n °t allowed due to the Pauli 



principle. 
The tl 

equations which, to leading order in a, explicitly are given by |^|, |l4], |l5| 



The three components, <j>^\. obey for large p the coupled angular Faddeev 



(fc& + ^ (a) ) ^»M a ) = 2a{-l) L cf p 2 Vi{p S ma) 



(4) 



2/ \ r-^C-^"^-^) 2 / • \ 2i fr\ 

«»(«) = -[ 2 P Viipsma) - v n \ , (5) 



where we defined v 2 = \ n + 4, the s-wave interactions Vi(xi) = Vjk{xi/ p,jk)jfi-, 
mp 2 k = mjmk/ (rrij + mfc) and the coupling terms 

cf = EE^1^t ' e = ' (6) 

where (p^ = arctan (mj. (r/ij + rrij + m/-) / rriimj) . The s-summation runs over 
the two possible intermediate spin couplings for the component j. The spin 
overlap matrix elements are diagonal for i — j, i.e. C" s , — 5 ss i and symmetric, 
le.CZ, = Cij s . 

The rescaled potentials p 2 Vi(p sin aj) approach for sufficiently large p the 
zero-range potentials, where the sensitivity to the shape disappears. Any po- 
tential with the same scattering length and effective range would then lead 
to results accurate to the order p~ 2 . We shall therefore for convenience use 
square well potentials |14 , where Vi^psmaC) = Vq^ (region I) when on < = 
&rcsiii(Rifj,jk/ p) <C 1 and zero otherwise (region II). 

The solutions to eq.(||) are then for a, > «q given by 
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L ^ (p,a i ) = AZ>P L (v,a i ), (7) 



P L (y, a) = cos L a ( ) sin z/ f a - ^) (8) 

\ oa cos a J I \ 2/1 

and for on < al to leading order in given by 



= 4° sin(^(0)a) - 2a(-l) i ^||cl: i) , (9) 



where the wave functions in in eq.(0) must be (fi^'l 
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Matching the solutions, eqs.(0) and (||), and their derivatives at on — off 
gives a linear set of equations for and . Large-distance physical solutions 
are obtained when the corresponding determinant is zero. This is the quanti- 
zation condition for A and the eigenvalue equation determining the asymptotic 
behavior of X(p). 

Input parameters. We use the neutron-neutron interaction which repro- 
duces the low-energy properties of the nucleon-nucleon system |jj . The neutron- 
core, n— 9 Li, effective interaction assumes that the spin of both 9 Li and n Li is 
|. We use the mean- field spin-orbit term 1 • s„ and we include a spin-spin 
term to differentiate between the two spin couplings for a given orbital angular 
momentum I. Such a spin splitting term is hard to avoid due to the strong 
spin dependence of the underlying basic interaction. Several phase equivalent 
parametrizations are possible. They differ in the number of two-body bound 
states occupied by the core neutrons which therefore subsequently must be ex- 
cluded in the computation || . The results are very close and we shall therefore 
here only use the shallow potentials without bound states. 

We use the radial form, exp(— (r/2.55fm) 2 ), for all neutron-core potentials. 
The strengths V s (lj n ) depend on orbital I and total neutron angular momentum 
j n for zero core spin, and for finite core spin also on the total spin s = 1,2. All 
possible s- and p-waves are included whereas other waves can be ignored to 
the accuracy we need. We use Vi(s 1/2 ) = -6.89 MeV, V 2 (s 1/2 ) = -7.51 MeV, 
Vi(pi/ 2 ) = -38.59 MeV, V 2 (pi/ 2 ) = -35.65 MeV, Vi(p 3/2 ) = 43.91 MeV and 
V2O53/2) = 46.85 MeV corresponding to virtual s-states at 0.247 MeV and 0.140 
MeV and p^-resonances at 0.75 MeV and 1.60 MeV. The p 3 / 2 -potential is 
purely repulsive thereby simulating the Pauli blocking by the core neutrons. 

The choice of these parameters is dictated by the accumulated information 
about the structure of 10 Li, i.e. a p-resonance at about 0.6 MeV, a low-lying vir- 
tual s-state and a small spin splitting of these states |^|, [TtJ . We also demand that 
the Li-binding energy is reproduced; we obtain 305 keV with the correspond- 
ing root mean square radius of 3.34 fm. The calculated fragment momentum 
distributions in break-up reactions also compare rather well with measured val- 
ues suggesting about 20% of the ^-configuration in the n Li ground state wave 
function §. The remaining part is in the s 2 -configuration. Limited ambiguity, 
although some, is then left for the parameters of these "realistic" interactions. 

The number of Jacobi polynomials in the basis expansion is carefully chosen 
to give accurate numerical results up to a distance, typically around 40 fm, 
where the asymptotic behavior is reached and from then on the asymptotic 
solutions eqs.(^)-(^) are used. The accurate low-energy continuum spectrum 
calculations require integration of the radial equations up to distances of the 
order of five times the sum of the scattering lengths. For the n+n+ 9 Li system 
this is about 200 fm. 
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Figure 1: The lowest angular eigenvalues A„ as functions of p for angular mo- 
mentum J* = | (dotted lines), and | + (solid lines) for n Li. The dot-dashed 
lines are the large distance asymptotic behavior. The levels marked ^ and 
| + also appear for those quantum numbers. The neutron-neutron interaction 
is from and the neutron- 9 Li interaction is given in the text. Maximum hy- 
perspherical quantum numbers up to K n 100 are used. 

Angular eigenvalues. The key quantities are the angular eigenvalues X n (p), 
which, divided by p 2 , are the adiabatic effective radial potentials. They depend 
on total angular momentum and parity as seen in Fig. 1, where we show the 
spectra for | (ground state) and | + (l~-excitation). The structure is com- 
plicated at small distances where avoided level crossings are seen. The lowest 
level has in both cases an attractive pocket, which is responsible for the ground 
state (| ) and several resonances. At large distance the structure is simpler 
as the hyperspherical spectrum is approached. In the computation we use the 
asymptotic behavior (shown in Fig. 1) of the three coupled s-waves, one from 
each Jacobi set. 

The 11 Li-states can be labeled as arising from couplings of the sx/2 and 
Pi/2 neutron-core states to angular momentum and parity 0^ and 1^, which 
in turn is coupled with the core spin 3/2 _ . The result for each parity 7r = ±1 
is three nearly degenerate states of 1/2*', 3/2*, 5/2* and one non-degenerate 
3/2*-state. This structure is reflected in the eigenvalue spectra where 3/2 ± 
contain additional levels on top of spectra similar to those of l/2 ± and 5/2 ± . 
The degeneracy would be complete without spin splitting of the neutron-core 
potentials. 
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Figure 2: The eigenphases corresponding to the lowest A„-values obtained after 
diagonalization of the S-matrix for J n = | (dotted lines) and | + (solid lines). 
The interactions are the same as in Fig. 1. 



Phase shifts, poles and the Efimov effect. We show in Fig. 2 the eigen- 
phases for the cases in Fig. 1 . For | we find rapid variations at smaller energy 
and one crossing of?r/2 at about 1.2 MeV. For \ and f " , none of them shown, 
we find rather smooth phase shifts indicating that the structure in § is due to 
the non-degenerate components. For | + four crossings of 7r/2 are seen in two 
pairs. In ± and | is found a similar structure, but one member of each pair 
is not present indicating the non-degenerate nature. 

For the lowest spins and parities we give in Table 1 the lowest S-matrix 
poles obtained by the complex energy method. We see the three times nearly 
degenerate l~-resonance at about 0.65 MeV with a width of about 0.35 MeV. 
We also see degenerate ± -resonances at 0.89 MeV with widths of 0.33 MeV 
and 0.43 MeV. In addition degenerate l ± -poles appear at 1.3 - 1.4 MeV, 1.6 - 
1.8 MeV and 1.9 - 2.1 MeV with widths of about 0.5 MeV, 0.6 MeV and 0.8 
MeV, respectively. Also poles of ± -type should be present for 3/2 T at these 
energies as seen by comparing with the model of zero core spin. They can not 
be distinguished numerically. 

The height of the effective radial barriers are 1.7 MeV, 0.9 MeV, 0.7 MeV 
and 0.6 MeV for ^ , | + , | + | + , respectively. The pocket and the barrier in 

this potential is absent for both i and | . These potentials are all attractive 
when the centrifugal barrier is removed. The lowest resonances appear around 
the barrier and their widths are consequently relatively large and rather sensi- 
tive to fine tuning of the interactions. Related cross sections would probably 
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Table 1: The real and imaginary values (E r ,T) (in MeV) of the lowest S- matrix 
poles E = E r — iT/2 for 11 Li for various spins and parities J 77 . The excitation 
energy E* = E r + 0.305 MeV. The interactions for the upper part of the table 
are the same as in Fig. 1. The lower part of the table are results for a model 
with s c = and the same average positions of the neutron-core resonances, i.e. 
energies of the Sin virtual state and the pi/ 2 resonance at 0.2 MeV and 1.22 
MeV, respectively. 





E r r 


E r r 


E r r 


E r r 


E r r 


1 - 

2 

3 _ 
2 

5 — 
2 


-0.305 


0.89 0.43 


1.37 0.51 
1.41 0.56 
1.36 0.49 


1.56 0.56 
1.60 0.61 
1.60 0.68 


1.98 0.65 
2.03 0.68 
2.01 0.72 


1 + 

1 + 

1 + 
2 


0.65 0.35 
0.68 0.33 
0.68 0.37 


0.88 0.33 


1.28 0.48 
1.33 0.50 
1.36 0.55 


1.74 0.64 
1.77 0.63 
1.74 0.64 


1.95 0.68 
2.08 0.71 
2.11 0.84 


0+ 

1+ 


-0.305 


1.00 0.37 


1.35 0.45 
1.40 0.59 


1.62 0.61 
1.59 0.63 


1.96 0.92 
2.02 0.81 


0" 

1- 


0.64 0.31 


0.92 0.39 


1.25 0.51 
1.46 0.53 


1.82 0.62 
1.76 0.59 


2.02 0.65 
2.08 0.67 



be relatively smooth. We have not attempted to reproduce a l~-resonance pre- 
cisely at 1 MeV II , but we still reproduce the properties almost within the 
experimental uncertainties. 

The relatively large number of additional poles could be due to the Efimov ef- 
fect, which is an anomaly in a three body system when the scattering lengths are 
much larger that the range of the interactions |j| . With increasing scattering 
lengths, the infinitely many poles of the three-body S-matrix move towards the 
point E=0. For very large but finite scattering lengths a number of poles must 
already appear close to zero. These poles originate from the long distance tail of 
the effective potential (oc Y2i=i 3 a i^JkP~ 3 ^ where a,i is the scattering length of 
the i-th subsystem) and they are not sensitive to the details of the interactions. 
Since there are no confining barriers for these poles, their corresponding widths 
must be rather large. 

In our case the Efimov condition is almost fulfilled, since the scattering 
lengths are much larger then the range of the interactions: a nn p,~^ + 2a cn \i~^ « 
50 fm. This must necessarily result in a number of broad resonances near the 
£■=0 point. 

Transition strengths and Coulomb cross sections. Dipole strength func- 
tions for different potentials are shown in Fig. 3 both for computed continuum 
wave functions and for plane waves. The interaction without spin splitting (ap- 



7 



2.5 




Energy (MeV) 

Figure 3: The dipole strength functions dB/dE are shown here for several 
potentials for comparison. The curves peaked at low energy are for the correct 
continuum wave functions and the broader curves are for plane waves. The 
interactions for s c — 3/2 are the same as in Fig. 1 and for s c — as described 
in table 1. JJH refers to the neutron-core potential from jl8|, where the result 
for the limited basis K max = 10 is from [Q. The measured curves are from 
[q 0-0- We normalize the data in |T(J to our sum rule value, while the absolute 
data from [g ^| is left unchanged. 

proximation with zero core spin) gives a distribution shifted about 100 keV 
towards lower energy compared to the result for the realistic full computation. 
A lower and broader peak is obtained for the potential from |l^| where the 
j» 2 -content of the three-body wave function is very small. For comparison we 
also show the result for the same potential with the limited basis employed in 
[ pjj . The low-lying l~-resonances enhance the strength functions at low ener- 
gies compared to the plane wave computation. The computed strength functions 
substantially exceed most of the data points (8[ ||, [l0| in the peak region around 
0.55 MeV. (Note that the data in [|[ ||] contain much less total strength.) A 
reduction could be achieved with higher resonance energy and larger width, 
but this would probably only be provided by a potential with much too small 
p 2 -content in the three-body wave function. On the other hand, the proper 
comparison with the different experimental results is not obvious. 

The Coulomb dissociation cross section is obtained by folding the virtual 
photon spectrum with the strength function. The precise expressions for the 
dipole approximation adopted here can be found in fis|| . The low-energy en- 
hancement places the strength where the virtual photon spectrum is relatively 
large and therefore necessarily also implies a large Coulomb cross section, since 
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the dipole is the dominating dissociation mode. The total cross sections for 
beam energies 42 MeV/A, 180 MeV/A, 280 MeV/A and 800 MeV/A on a lead 
target are 4.487 b, 2.128 b, 1.429 b and 0.971 b, respectively. This is in agree- 
ment with the available measured values ||, but substantially larger than 2.586 
b, 1.166 b, 0.933 b and 0.657 b found by using plane waves. For 280 MeV/A 
we find 1458 b, 1501 b, 1283 b, respectively for the potentials without spin 
splitting, with small p 2 -content with complete and limited basis. 

Conclusions. A recently formulated method was used to compute low-energy 
three-body continuum spectra for n Li for various angular momentum states. 
The angular part of the Faddeev equations are treated numerically at short dis- 
tances, whereas the large distance behavior of eigenvalues and eigenfunctions is 
computed essentially analytically. Combining the results from these two regions 
allow accurate computations at large distances. 

We employ interactions that reproduce reasonably well the low energy con- 
tinuum properties of the two-body subsystems and the fragment momentum 
distributions in break-up reactions. The correct spins of 9 Li and 11 Li are used 
implying level splitting and fine structure beyond the approximation of zero 
core spin. A l _ -resonance at about 0.65 MeV with a width of 0.3 MeV is found 
in fair agreement with a recent measurement. Also several other s-matrix poles 
were found and related to the Efimov effect. Strength functions, Coulomb disso- 
ciation cross sections, phase shifts and S-matrix poles for J"" = ^ , ^ , ^ are 
computed. A proper comparison of measured and computed strength functions 
would be rewarding either by showing the limitations of, or by selecting, the 
most correct three-body model. 
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